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ABSTRACT 


Errors  in  the  division  of  power  by  a variable-power  divider  using  phase 
shifters  and  a pair  of  3-dB  couplers  are  determined  as  a function  of  an  incor- 
rect setting  of  the  phase  shifters.  The  insertion  phase,  of  the  phase  shift- 
ers, is  assumed  to  have  a Gaussian  distribution  about  the  desired  value  and 
the  resultant  probability  distribution  of  the  power  division  is  calculated. 
Curves  are  given  for  the  probability  that  falls  within  a given  error. 

These  curves  are  computed  from  the  probability  density  function  of  P^  which 
is  also  shown  graphically. 
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Error  Analysis  of  a Variable-Power  Divider 
A variable-power  dividing  four-port  junction  consisting  of  a magic  "tee”, 
a 3-dB-short-slot  coupler,  and  two  latching-ferrite  phase-shifters,  as  indi- 
cated schematically  in  Fig.  1,  has  output  signals  given  by  the  following 
equations: 


where  <j>^  and  ij>^  are  the  insertion  phase  shift  of  the  latching-phase  shifters. 
These  equations  assume  the  magic  tee  and  3-dB  coupler  are  perfect  lossless 
devices,  an  assumption  that  is  reasonable  in  view  of  the  comparative  error  in 
establishing  and  (f^. 
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Fig.  1.  Variable-Power  Divider. 
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Equations  (1)  and  (2)  can  be  put  in  a somewhat  simpler  form  if  we  assume  $2 
varies  from  ~ to  77  and  <f>,  varies  from  0 to  that  is  <p^  has  a minimum  inser- 

TT 

tion  phase  of  ~ radians.  In  addition  we  let  £ ■ ~4)^  - tt/2,  and  consider  only 

the  power  delivered  to  the  output  ports  since  the  output  signals  are  always 
in  phase  (because  we  also  assume  the  same  insertion  loss  for  each  phase 
shifter)  and  our  interest  centers  on  the  device  as  a power  divider.  Hence 

PA  - sin2  | (3) 

and 

PB  *»  cos2  | (A) 

The  length,  T,  of  a pulse  applied  to  the  ferrite  phase  shifter  determines 
its  Insertion  phase.  The  relationship  between  insertion  phase  and  pulse 
width  is  non-linear;  it  varies  as  a function  of  temperature  and  aging,  and 
each  device  is  also  somewhat  different  from  all  other  "identical"  devices.  The 
problem  addressed  in  this  note  concerns  determining  the  error  in  P.  and 
when  a polynomial  representation  such  as; 


N 


is  used  to  set  the  phase  shifters,  in  any  variable  power  divider  (VPD) , to 
obtain  a desired  division  of  power. 

Toward  this  end  let  us  first  determine  a method  of  obtaining  (5),  or  its 
inverse 


2 


(5a) 


n-0 


and  what  probability  density  function  (pdf)  we  should  assign  to  the  dependent 
variable.  Although  it  is  customary  to  determine  T for  a given  power  division 
by  using  (5)  and  (3),  T is  ultimately  applied  to  the  control  equipment  so  as 
to  obtain  the  desired  power  division.  Hence  it  is  more  practical  to  consider 
T as  the  independent  variable  and  use  (5a)  and  (3)  or  (A)  to  evaluate  the 
error  in  PA  and  Pfi. 

In  particular,  let  us  determine  the  insertion  phase,  d>  of  a large  num- 

mn 

ber  of  phase  shifters  say  M«  AO  or  more,  (m  * phase  shifter  number,  n ■ drive  pulse 
duration)  for  a f in..te  set  of  error  free  applied  control  pulses  • .xN)  and  sev- 

eral temperatures  and  frequencies  in  the  operating  range.  Assuming  $mn,  for  a partic- 
ular value  of  n * i,  to  have  a Gaussian  distribution  about  the  mean  value,  4> . • 

_ M __  1 

(cf>.  ■ (1/M.)  E d>  . ) , the  standard  deviation  a.  and  $.  are  computed.  (In  a 
I _ , ml  l i 

m*l 

most  exact  analysis  one  would  attempt  to  determine  the  actual  frequency  func- 
tion, or  pdf,  but  the  assumed  Gaussian  and  perhaps  a uniform  distribution 
should  establish  an  interesting  bound  as  will  be  discussed  later.)  The 
values  of  ^ and  are  then  fit  to  a pth  degree  polynomial  by  a suitable 
regression  analysis  and  the  an  of  (5a)  are  determined.  It  may  also  be  useful 
to  determine  the  b , in 

M 

0 “ V™  (6) 

m-0 
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in  the  same  manner.  Hence  (5a)  and  (6N,  together  with  assumed  Gaussian  dis- 
tribution, define  the  statistical  performance  of  the  phase  shifters  as  a group. 
For  a Gaussian  distribution  of  <J>(t^),,  for  fixed  , the  linear  combin- 


ation 


5 - 4>(t4)  - 4>(xm)  (7) 

is  also  a Gaussian  distribution  with  mean  £ = <p„  - i and  variance 

it  m 


It  remains  to  find  the  pdf  of  p^  (u)  » and  then  the  probability  that  PA 
lies  between  P^+  A PA  and  P^-AP^.  Toward  this  end  we  recognize  that  the  prob- 
ability density  function  of  £ is  given  by 


(x)  = 


. I(2^)2 

2 ' o ' 
/2tt  a 


(8) 


The  relationship  between  P^,  or  Pg  and  £ is  obtained  from  (3)  and  (4)  and 
rewritten  here  as 


P„  - 


11- 
2 " 2 cos^ 


1,1  - 

2+2  cosC 


1-P. 


(9a) 

(9b) 


The  probability  x^  < P^  < X£  is  given  by 

x2 


Pl*l  < PA  < 


X21  -J 


PA(u)du 


(10) 


where  x-j_  and  X£  are  assumed  values  of  PA< 
Sketching  (8)  and  (9b)  (see  Fig.  2a  and  2b), 


Aik 


Fig.  2a.  Probability  density  function  of  £. 


Fig.  2b.  Distribution  of  P 


it  is  apparent  that  PA  is  a periodic  function  of  £ and  obtains  a value  less 

than  x for  all  £ between  -£j  + 2n'rr  and  £ + 2mt  where  n is  any  integer 

X X 

between  -»  and  °°.  Hence  from  Figs.  2a  and  2b  we  can  write 


- cos-1[l-2x]  . (14) 

If  we  limit  K to  the  range  0 £ £ tt  all  interesting  values  of  PA  will  be 

obtained;  i.e.,  0 <_  PA  <_  1.  We  should  also  recognize  that  a in  (8)  is  »0.1. 

Consequently  for  0 < X < ff/2  (0  < x < 1/2),  only  the  n *=  0 term  in  (13)  will 

be  significant.  For  ir/2  < E,  < tt  (1/2  < x < 1),  the  n ■ 1 term  must  also  be 

taken  into  account.  Specifically 

Z- 

P[PA  < x]  **  r*  p^(u)du  (13a) 
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for  0 < x < and 


} HI 

a 


2-nr 

[PA  < *'  * J p^(u)du  ^ p^(u)du 


(13b) 


for  < x < 1.  Note  that  (13b)  neglects  the  contribution  for  -£j  < £ < 0 

and  2h  < ^ < 2tt  + C because  they  are  negligible. 

X 

Recalling  the  definition 


pA(x)  * Lim 


P[x  < PA  < x + AP  1 


AP.-K) 

A 


(p[pA  < X]) 


and  using  (13a) , we  have 


V'X)  -W  1?  “ P^x>  “~T x”  • 


(16a) 


Using  (13b)  we  have 


d<-v  1 


Vx)  -Pf,(y  dT  -pe<2*-V  -sr-  • 2 * * i 1 


(16b) 


From  (14),  we  have 


^Sc  „ __1 

dx  /x(l-x) 


Substituting  (17)  and  (8)  into  (16a)  and  (16b)  we  consider  the  combination 
resulting  from  (16a)  and  (16b)  to  be 
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■w  a—— aaaawaa  n nil  mi  SjSMMBBf 


(18) 


It  is  clear  that  only  the  first  term  on  the  right  hand  side  will  contrib- 
ute for  0 < £ < tt.  The  second  term  contributes  when  £^  and  £ approach  0 and 
the  third  term  contributes  only  when  £ and  £^  approach  tt.  As  we  would  expect 

PA(x)  has  its  maximum  value  for  £ « X-  The  probability  density  function 

2— 

PA(x)  is  sketched  in  Fig.  3 with  = sin  £ as  a parameter. 


Ill  - 6-  11647 1 


Note  that  the  pdf  has  singularities  at  x «*  0 and  1.  However,  it  is  an 
integrable  singularity  and  the  area  under  curves  shown,  over  the  range 
0 < x < 1,  is  finite  and  equals  1 as  it  should.. 


7 


For  example  let  us  evaluate  the  area  over  a region  (0,e)  about  the  x - 0 sing 


ularity,  i.e.. 


f„ 


p[?A  < e ] « / PA(x)dx 

0 

where  £ ■*•  0.  Using  (18)  we  have 


(19) 


P[PA  < e]  « 


a/Ii t 


-|(^)2+  -|(^)2ld 

• ±__e j dx  (20) 

/x 


where  we  have  set  ^ = 2/x  (see  (14))  and  /x(l-x)  - v^T  Rearranging  (20) 


P[PA  < e] » 


i/?f  . -2!s*&  -2±lM 

2\aJ  r a 2 


gVZtT 


l 


4*  e 


& 


dx.  (21) 


For  x 0 the  numerator  of  the  exponent  can  be  approximated  as  indicated  in 

(22) 


?,2 


P[PA  < e] 


“i(f) 


cm 


£ 

/ 


„ 4x  , 2x2  , 2x£2 

7 J 

o a a 


dx 


(22) 
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-Kl  i 


o/ffr 


1/2  4(2-g)e3/2  4e5/2 

' 3a2  5o2 


] 


(23) 


Now  if  e < < a (recall  0 < £ < tt)  , the  second  and  third  terms  (23)  can  be 
neglected  and 
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liMMiiaaaittiiiM 


llnl  iv'l 


P[PA  < e]« 


-Ki)2 

o/Jrr 


Notice  that  (24)  gives  a finite  probability  indicating  that  integration  over 
the  singularity  at  x ■ 0 does  indeed  give  a finite  value.  With  c m 0.1  and 

X - 0.2 


P[PA  < £]  « 1ST  (25) 

where  e < .001.  Conversely  P[PA  > e]  ■ 1 - 1.5/eT.  This  result  will  be 
checked  in  the  next  section, 
milative  Distribution  of  P, 

A 

Using  the  foregoing  expression  (18),  and  using  numerical  integration. 


A 

< x]  ^ 


PA  (u  ) du 


was  calculated  for  o * 0.05,  0.1  and  0.2  and  PA  “ 0.01,  0.5,  0.99.  PA  is  the 
desired  value  of  PA;  it  is  not  necessarily  the  mean,  or  average,  value  of  PA. 
The  results  are  shown  in  Figs.  4,  5 and  6.  The  P[PA  - x <■  PA  < PA  + x]  is 
plotted  versus  magnitude  of  Ax  ■ (P^  - x)  in  Figs.  7,  8 and  9. 

Notice  that  for  PA  “ .01  (Fig.  4),  which  corresponds  to  $ * .2, 

P[PA  < .001]  s .07  which  agrees  reasonably  well  with  the  approximate  relation- 
ship given  by  Eq.  (25).  It  is  also  important  to  note  the  asymmetrical  character  of 
the  curves  with  respect  to  x’  ■ PA  ■ .01  the  desired  value  of  PA<  For  a « .2 
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P 0.01  - Ax  < P<  0.01+Ax 


«V+66  0 > d > *V  -66  O I d 


r *vto5*?K*1 


**yfr,  .<.^1  *.'.■  tf*-  < »} 


the  mean  value  of  P^,  P^,  is  » 0.011! 

The  knee  in  the  curves  shown  in  Fig.  7 and  Fig.  9 i3  due  to  the  asymmetry 
in  the  cumulative  distribution.  In  particular  the  difference  in  slope  for 
x > P^  to  that  for  x less  than  P^  causes  the  abrupt  change  in  the  slope  of  the 
curves  (Figs.  7 and  8)  since  they  are  computed  from  the  cumulative  distribu- 
tion. 

Uniform-1 , Distributed  Errors 

Let  us  next  assume  that  the  errors  in  Eq.  (5)  or  (5a)  have  a uniform 

distribution  instead  of  the  previously  assumed  Gaussian  distribution.  It  can 

be  shown  that,  if  4>(T^)  is  distributed  uniformly  within  ^(t^)  + <f>0* 

£ * 4>(t0)  - 4> (T  ) will  have  a mean  value  £ ■ 4>(t0)  “ <Kt  ) and  it  will 
X/  in  Xj  in 

have  a triangular  distribution  between  X + 2<bf  . In  particular,  the  pdf  be- 
comes 

pF  ■ <1+iS)/2fo!  -24>„  < “ - ? < 0 

’ o 

(28) 


■ <2  - if)''2*,1  0 < “ - s < 2*„ 

To 

Using  (13)  and  (14) 

P[P,  < x]  - — 

4*o 

when  £ < X (29a) 

K;  - Z - 2d»0,  ex  “ o if  K > 2(|»o.  (29b) 

If  I < 2<(,o, 


(2(f>o“^+u)du 


x 

+ ‘1/ 


(2<t>o-f-u)du 


(29) 


16 


-.  .-t^t^v^  j-.'wr  r,-*^?^'w^17^j^'iV7^^>- , 


s^w*rw?.y  wffiirr'' 


5 ' - 0 
^x 

Cl  - 1 if  CX  < 2<j>o  - T 

£l  « o if  Cx  > 24,o  - T 


(29c) 

(29d) 

(29e) 


when  5 > 5 

x — 


P[PA  < X] 


t?f 


(2<j»0-T-u)du  + P[PA  < sin2  \\ 


(30) 


_ 2 
where  L,  < ir  ~ 2<j>o  which  corresponds  to  x < cos  $g.  It  is  also  necessary  that 

5 <5+2 6 . In  order  to  calculate  P[P.  < x]  when  ■tt— 2<J>  < 5 < it,  we  can 
^x  — ro  A o 


use  (29)  for  5X  < 5 and  an  appropriate  modification  of  (30).  Because  this  is 
not  a particularly  interesting  case  (i.e.,  5 > n " 2^)  it  will  not  be 
treated  here. 

Integrating  (29)  and  (30)  requires  attention  to  the  conditions  indicated. 

Let  us  first  consider  those  cases  where  24  < 5 < fr  - 2<j>  . For  E < 5,  let 

o ox 

c’x  mT  - 24»0;  then 


P[PA  < x]  - 

For  5x  > I 

P[PA  < x]  - 

When  5 < 2<j>  ,we  have 
o 

P[PA<  x]  - 


8(j) 


7<5x-5’x)2. 


i + 7T2  ' 5*  + » «,-«• 

84 


K\ 


X X 

2d  2 


(31) 


(31a) 


(32a) 
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for  0 < £ < £ and  £ < 2 <b  -£.  There  are  two  possible  conditions  to  consider; 

x x o 

I . e . , X < 2<j)o~£  and  £ > 24>q-£.  Taking  the  latter  case  first 


r‘PA  < *1  --,72  + “4  + C-O 

2<^o  8*0 


for  2<}!0-£  < ^ < £ and 


(32b) 


p[PA<x1-i  +-72  WxKx-K\->  * KxWi 

84> 


(32c) 


for  £ < £ < £x  , where  £x  " 2<J>o  + £. 

Now  considering  the  case  £ < £ ' , we  have 


P[PA  < X]  - -L  M05x  - 5X2  + E2] 


(32d) 


for  £ < £x  < £'  and 


P[PA  < x] 


-4  I?’*2  - 2E2  * 25;ex-ex2]. 


(32e) 


when  £'x  < £x  < £"x  . 


A plot  of  P[PA  < x]  for  £ - tt/2  (P'a  » 0.5)  and  4>q  - 0.1  radian  is 
shown  in  Fig.  10.  Notice  the  similarity  between  this  curve  and  that  for 
a a .1  in  Fig.  4.  As  should  be  expected,  uniformly  distributed  errors  with 
a maximum  value  4>  “0,  the  RMS  value  of  errors  with  a Gaussian  distribution, 
are  more  likely  to  produce  the  desired  value  of  PA»  This  present  analysis 
shows  the  possible  effects  should  the  actual  errors  be  distributed  either  as 
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Gaussian  or  a finite  uniform  distribution.  Using  the  foregoing  analysis,  it 
is  also  interesting  to  note  that  for  <j>Q  - 0.1  radian  and  uniformly  distrib- 
uted errors,  an  attempt  to  set  < 0.01  will  be  successful  only  50%  of  the 
time. 

Conclusions 

Should  the  measured  data  show  that  the  errors  have  a Gaussian  distribu- 
tion and  they  have  an  RMS  error  ■ 5°  (*  .1  radian)  one  can  expect,  from 
Fig.  8 (P^  ■ .99),  that  80%  of  the  time  the  power  delivered  to  the  desired 
port  (say  port  A)  will  be  within  0.05  dB  of  the  maximum  power  available.  At 
the  same  time  the  power  delivered  to  Pg  will  be  more  than  17  dB  below  power 
into  the  variable  power  divider.  Similarly  from  Fig.  7,  for  equal  power 
division  (P  A "0.5),  80%  of  the  time  the  P^  will  be  within  0.3  dB  of  the 
desired  value.  From  these  examples  we  see  that  only  when  trying  to  minimize 
P^  (or  Pfi)  does  the  variation  of  <|>(t)  effect  the  power  division  significantly. 
Use  of  Eq.  (32a),  with  £ “ 0 (P'A  = 0),  indicates  that,  with  the  assumed  5° 

RMS  error  in  4>,  P^  will  be  less  than  -30  dF>  (referred  to  the  input  power) 
with  a probability  of  0.62.  This  operational  performance  may  have  question- 
able value  in  a nulling  nntenna  when  null  depth  greater  than  20  dB  is  desired. 
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